function Costs = generateCosts( nbOfCities )

C =nbOfCities;
sigma = .5;
rho = .05;
SigmaInd = sigma.^2 .* [1 rho; rho 1];
Clients = mvnrnd([0.5 0.5], SigmaInd, C);
Clients = (Clients +2).*50;

XDist = zeros(C, C);
for i=1:C
    for j=1:C
        XDist(j,i) = sqrt((Clients(i,1) - Clients(j,1))^2 + ...
            (Clients(i,2) - Clients(j,2))^2);  
    end
end
Costs = XDist;
end

